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O ■ The probability of a projectile nucleon to traverse a target nucleus without inter- 

action is calculated for central Si-Pb collisions and compared to the data of E814. 
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O ! The calculations are performed in two independent ways, via Glauber theory and 

O ■ using the transport code UrQMD. For central collisions Glauber predictions are 

about 30 to 50% higher than experiment, while the output of UrQMD does not 
show the experimental peak of beam rapidity particles. 
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^ \ The dynamics of a heavy ion collision is a complicated many-body problem. 

It is the task of appropriately designed experiments to isolate one particu- 
lar aspect of the dynamics and elucidate its physics. Wounded nucleons are 
one of the open problems. In a heavy ion collision a nucleon may undergo 
a sequence of collisions, w^hich follow so rapidly one after another, that the 
nucleon is no more in its ground state, even not necessarily in any definite 
excited baryonic resonance (like A or N*). We will speak of a wounded nu- 
cleon. In a next encounter with another nucleon this wounded nucleon will 
not interact with the free space NN cross section af^^ but with an effective 
one (Jeff. Is it possible to determine cTefr from experiment? 
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The E814 collaboration has designed an experiment which may answer this 
question [1]. At the energy of 14.6 GeV/nucleon, the projectile ^^Si collides 
with Al, Cu and Pb and the beam rapidity nucleons are studied as a func- 
tion of the centrality of the reaction (controlled by a measurement of the 
transverse energy). The beam-rapidity nucleons belong to the projectile and 
have not lost any energy in the reaction. Of course, in a peripheral reaction 
one always has the so called spectator nucleons, which never hit the tar- 
get nucleus. They are not interesting for our purpose. However, in a central 
event, e.g., for Si on Pb, one still sees beam rapidity nucleons, i.e., nucle- 
ons of Si which go through the target nucleus without energy loss. The 
target nucleus is transparent for these projectile nucleons. We expect the 
transparency of a heavy nucleus like Pb to be small. Indeed, the "survival 
probability" S for a projectile nucleon to pass through the target without 
any inelastic interaction has been measured to Sexp — 3.5 x 10~^ for central 
Si-Pb collisions [1]. Does this result contain information about wounded 
nucleons and effective NN cross sections? 

We proceed in the following ways: The transparency is calculated without 
assuming any exotic phenomena and then the quality of agreement with 
experiment is judged. We use two complementary methods of calculation 

• Glauber theory, in which the physics and formalism is simple and trans- 
parent, though not always realistic; 

• Cascade calculation (UrQMD), which represents present-day state-of-the- 
art and may be compared to earlier calculations (Bass et al. [2]). 

The survival probability S for a projectile nucleon to pass through the target 
without inelastic interaction is calculated in a Glauber-type approximation 
(straight trajectories, frozen nucleons) as 

Saib, (Jeff) = / d^sTpib- s)e-''<'^^*^*^^, (1) 

where 6 is the impact parameter of the nucleus-nucleus collision, Tp{b) and 
Tt{b) are the thickness functions {T{b) — J dzp{b, z), J d?xp{x) = 1) of projec- 
tile and target, respectively. The straight line geometry is certainly justified 
for the through-going nucleons. The use of "frozen" nucleons and the ne- 
glect of any other degrees of freedom, like mesons, need justification. We 
consider a target nucleon and estimate the time At in its rest system which 
has elapsed between the arrival times of the first and the last projectile nu- 
cleons: At < 2i?p/7p where 2Rp ~ 7 fm is the diameter of the Si projectile 
and 7p ~ 16 is the Lorenz factor for the experiment under consideration. 
The time At ^ 0.5 fm/c is rather short: (i) The target nucleon has not 
moved significantly in space (frozen approximation is good), (ii) Accord- 
ing to the uncertainty principle the intrinsic excitation energy is uncertain 
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Fig. 1. Survival probabilities SG{b,cTcs) of a beam rapidity nucleon calculated in 



Glauber approximation for a^s = c". 
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30 mb and cJefr = a, 



NN 
tot 



39 mb (solid 



lines), and the experimental value Sexp for the survival probability for = 15.5 
GeV. The correlation function P{Et, b) between transverse energy Et and impact 
parameter b is (i) calculated in Glauber approximation, PciEt, h) (continuous line), 
and (ii) deduced from the results of UrQMD, Pu{Et, h) (histogram), both plotted 
as bP in order to represent da/db. 

by > (At)"^ ^0.4 GeV and does not allow the assignment of a defi- 
nite excited state, (iii) Secondary hadrons are not yet formed since typical 
formation times are of order 1 fm/c. 

The value of the transparency S in Eq. (1) depends on the effective cross 
section cXefF. Using Saxon- Woods parameterizations for the densities pt and 
pp with the surface thickness a = 0.52 fm for all nuclei and the half-density 
radius ta such that the root-mean-square radius of the nucleus equals the 
charge radius [3], the survival probability S{h, a^s) is calculated for two val- 
ues afj^ = 30 and a'^J^ = 39 mb which correspond to the inelastic NN cross 
section and the total one at this energy, respectively. The results are dis- 
played in Fig. 1. The experimental point, measured at transverse energy 
— 15.5 GeV is also shown in the figure. It is obtained from the measured 
mean multiplicity (Mc) of beam rapidity protons by 



5e.p(^t) = {M,){Et)/Z, 



(2) 



where Z = 14 is the number of protons in Si. We have assumed - as the 
authors of the experiment do - that the highest ii^^^^-bin corresponds to a 
central collision which is assigned a value b — 0. Then the experimental 
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value is close to the curve = 30 mb. In fact the equation Saib = 0, aes) — 
Scj.p(E^ = 15.5 GeV) leads to a value a^s = 31.1 ± 0.7 mb. This result repro- 
duces a similar calculation using uniform density distributions by the E814 
collaboration [1], who have concluded that a^s = (^^^ within error bars and 
no anomaly being visible. 

The crucial step in the argument is the assignment of 6 = to the highest 
Et-hin. Actually, a given value of Et selects a distribution of impact param- 
eters only within a band Ab{Et) which is quite large. We study the relation 
between Et and b and Ab in the form of a probability distributionP(£'t, b), 
which gives the probability that values of 6 contribute to events with a given 
value of Et. We normalize it as / dEtP{Et, h) = 1. With this function and the 
differential inelastic heavy ion cross section d^ain/cPb, one can obtain the 
dependence of the survival function S{Et, cxefr) as a function of the trans- 
verse energy 

, Jd%Sib.a,s)P{Et,b)daiJd% 
b[Et,aes)^ — , (J) 

dain/dEt 



where 

^./,^JP(^„,)^. (4) 

We will use Eq. (4) to determine P{Et, b) from a comparison with the mea- 
sured transverse energy distribution duin/dEf With a reasonable Gaussian 
parameterization of PciEt, b) for the Glauber calculation (or with Pu{Et, b) 
from a UrQMD calculation, see below), S{Et, (Tcff) can be determined from 
Eqs. 1 and 3 without ambiguity. We discuss our parameterizations. 

The inelastic cross section <Pain/d% is taken from the folding model 

dar 



d% 



(6) = 1 - expi-aZ^'ApAt J dhT,{b- s)Tt{s)]. (5) 



The parameterization of P{Et, b) is more model dependent, and we cal- 
culate it in two ways: (i) in the formalism of Glauber theory, and call it 
PciEt, b), and (ii) deduce it from results of the cascade calculation, Pu{Et, b). 
To calculate Pq we assume a Gaussian parameterization 



which satisfies the normalization condition. We make the usual assump- 
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Fig. 2. Transverse energy distributions in the rapidity range —0.5 < r? < 0.8 for 
Si+Pb collisions at 14.6 GeV/nucleon. Data points (dots) are taken from [7]. The 
solid line is calculated within the collision model with parameters eo, a; determined 
by a fit to the Pb data. The UrQMD calculation (histogram) contains a simulation 
of the leakage of the TCAL counter as explained in the text. 

tions [4,5] for the functions Et{b) and at{b) 

Et{b) = 7V(6)eo, (7) 
al{h) = N{h)eluj. (8) 

Here N{h)\s calculated in the "collision model" [4] 

Nib) = a^^^'AA, J d%T,{s)T,{b- s), (9) 

which equals the mean number of NN collisions in a projectile-target in- 
teraction with impact parameter b. The proportionality constant eo between 
the observed transverse energy Et and N{b) depends on the dynamics of 
hadron production, but also on the experimental set up (chosen rapidity 
interval and acceptance of counter). It will be a fit parameter. The ansatz 
oc N{b) corresponds to the hypothesis that the fluctuations in N{b) are of 
statistical origin. It is known, however, that the proportionality constant uj 
in Eq. (8) depends strongly on the rapidity interval of the accepted particles 
[4]. The reason is not clear [6]. We take a; as a free parameter. 

Using expressions (5) to (8) we have calculated dain/dEt and have varied 
the parameters eo and ou until the data for Si-Pb are fitted. The result is 
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shown in Fig. 2 with a,^^ = 30 mb, eo = 0.067 GeV, uj = 7.8. A calcula- 
tion of N{b) in Eqs. (7-8) within the wounded-nucleon model [5] gives a 
similarly good fit to the data for dain/dEt. 

According to Eq. (3), P{Et, h) determines the integration region in impact 
parameter h, which contributes to the integral for a given value of Et. Fig. 1 
shows the distribution hPc{Et,h) for a central collision {E'j: = 15.5 GeV) 
derived in the Glauber formalism, and the corresponding one, bPu{E^, h), 
calculated (see below) with the help of the code UrQMD [8]. Both distribu- 
tions agree within statistics, with a peak around 1.7 fm with a half width of 
about 2.5 fm. Using the functions Scih, (Tcfr) and Pg{E^, h) as shown in Fig. 1, 
the calculated value of Sg{E^, aes) gets contributions from a considerable 
range of values of b. We find 

S{E^, a,s = SOmh)/ Sexp = 3.47 (10) 

The calculated transparency is too large by a factor of three.. Should one 
rather use = cr^^ instead? We think yes: Although elastic scattering 
changes the rapidity of the projectile nucleon only marginally, it increases 
the transverse momentum of a nucleon, and an elastically scattered nucleon 
is not accepted within the experimental window pt < 0.3 GeV with a large 
probability. The results of the Glauber-type calculations are compared with 
experiment in Fig. 4. We discuss this figure after we have described the 
second approach. 

The UrQMD calculations presented in this paper are based on the standard 
version 1.1 [8]. As mentioned in [1], due to leakage in the TCAL [7], there is 
less than half of the Et visible. This has two consequences: 

(i) The Et scale of the data cannot be directly compared with that of the 
cascade calculation, (ii) The data have additional fluctuations, visible as a 
flatter slope of the E814 i?t-distribution beyond the knee. Simulating the 
leakage in the Monte Carlo by simply selecting only 37% of the particles 
hitting the TCAL fixes both problems. The resulting ^^^-distribution agrees 
quite well with the data (Fig. 2). In order to compare to the E'j-selected pro- 
ton rapidity distributions of E814, we have used the absolute cross section 
given there - 593, 102 and 6.9 mb - for the three centrality bins shown in 
Fig. 3. The mean impact parameters for these bins are 3.6, 2.1, and 1.6 fm, 
resp., in UrQMD. 

Fig. 3 shows the rapidity distributions for protons + neutrons for the three 
Et bins. The data points are shown by the dots, while the two histograms 
represent the results of the UrQMD calculation with an upper cut in pt 
at 0.3 GeV/c. The two histograms show results assuming either perfect 
acceptance or taking the E814 experimental acceptance cut into account 
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Fig. 3. The rapidity distributions da/dy for protons and neutrons with pt < 0.3 
GeV/ c for three intervals of transverse energy Ef. The data are indicated by dots, 
the calculations using UrQMD are represented by histograms. In each figure the 
histogram with higher values corresponds to ideal acceptance, the lower, hatched 
histogram is obtained with E814 acceptance cuts. Beam rapidity is at y = 3.44, 
where the data show a clear peak while the calculation does not show one for the 
two higher Et bins. 

(hatched area). The two histograms are rather similar in the region of beam 
rapidity {ys = 3.44). While calculations agree with the data in the interval 
6 < £'( < 11 GeV, there is a clear discrepancy in the region around beam 
rapidity for the two other Et intervals (middle and lower picture): The cal- 
culations do not show any peak, while the data do. 

Fig. 4 shows a comparison between experiment and our two calculations for 
the multiplicity of beam rapidity protons. In the analysis of the data and in 
the output of UrQMD, two different criteria are employed to characterize 
"beam rapidity protons": (i) all protons with y > 3 (black dots for the ex- 
periment and triangles for the UrQMD calculation); (ii) twice the number of 
protons with ys > 3.44 - assuming the distribution to be symmetric around 
ys (open circles for the experiment and black squares for the UrQMD cal- 
culation). While the two criteria do not give very different results for the 
experiment, they lead to very different values for UrQMD up to a factor of 
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Fig. 4. Comparison of experimental and theoretical results for the multiplicity of 
beam rapidity protons as a function of transverse energy Et for the reaction Si+Pb 
at 14.6 GeV/ nucleon. Two methods are employed to characterize "beam rapidity 
protons" in experiment and in the output of UrQMD: (i) protons with rapidity 
y > 3 (black circles for the data and triangles for UrQMD) and (ii) multiplicity of 
protons with y > 2/^(3.44) multiplied by a factor 2 (open circles for the data and 
squares for UrQMD). The dashed and solid lines represent Glauber calculations 
with (Teff = (Jin = 30 mb and cTeff = <^tot = 39 mb, respectively. 

five smaller for the second method. This reflects the fact already observed 
in Fig. 3 that UrQMD does not produce a peak at beam rapidity for the 
more central values of Et. It is therefore not meaningful to compare data 
and results from UrQMD for Et > 10 GeV. For values oi Et < 10 GeV cal- 
culation and data agree to better than a factor of two - the calculations are 
above the data for Et — 5 GeV and below the data at Et — 10 GeV. The re- 
sult of the Glauber calculations (solid curve for aes — (^tot — 39 mb, dashed 
curve (Tcff = cTin = 30 mb) are nearly always above the data. The solid curve 
agrees within 30-50% for central values of Et, but agreement worsens to- 
wards smaller values of Et. At Et = 5 GeV the Glauber calculation is a factor 
3 higher than the data. Note that in this region also the results of UrQMD 
lie above the data by a factor of 1.5. These discrepancies may partly have 
the following origin: Small values of Et correspond to peripheral reactions 
with several nucleons being spectator particles. All spectator nucleons are 
included in the calculations; spectator nucleons bound in deuterons, alpha 
particles etc., however, are not contained in the data. 

In summary, the two theoretical approaches, Glauber theory and UrQMD, 
applied to the transverse energy dependence of the nuclear transparency, 
have yielded only a semi-quantitative understanding of the data. At small 
values of Et (peripheral reactions) calculations are systematically above the 
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data. In the interesting regime of central collisions, Glauber theory predicts 
values of 30 to 50% above the data, while UrQMD shows clearly no peak of 
beam-rapidity particles so that the results are only an upper limit. 

Although the theoretical situation is still far from being satisfactory the 
analysis of beam-rapidity particles in heavy ion collisions may potentially 
be a good way to investigate the question of effective cross sections of 
wounded nucleons. 
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